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We show a new physical phenomenon expected for the ratio R p a of the unintegrated gluon dis- 
tribution of a nucleus over the unintegrated gluon distribution of a proton scaled up by the atomic 
factor A 1 / 3 in the fluctuation-dominated (diffusive scaling) region at high energy. We calculate the 
dependence of R v a on the atomic number A, the rapidity Y and the transverse gluon momentum 
k±. We find that R p a exhibits an increasing gluon shadowing with growing rapidity, approaching 
l/A 1 / 3 a t asymptotic rapidities which means total gluon shadowing, due to the effect of gluon num- 
ber fluctuations or Pomeron loops. The increase of R v a with rising gluon momentum decreases as 
the rapidity grows. In contrast, in the geometric scaling region where the effect of fluctuations is 
negligible, the ratio R p a shows only partial gluon shadowing in the fixed-coupling case, basically 
independent on the rapidity and the gluon momentum. 
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There has been a tremendous progress towards understanding gluon number fluctuation or Pomeron loop effects 
in the high energy evolution in QCD over the last two years. QCD evolution equations have been established, 
sometimes called "Pomeron loop equations" P, H, Q, which take into account fluctuations and a relation between 
high density QCD and reaction-diffusion processes in statistical physics has been found which allows us to obtain 
universal, analytical results for scattering amplitudes in the fluctuation-dominated regime at high energy 0, 0, H[ . 
The main result , as a consequence of fluctuations, is the emerge of a new scaling behaviour for the dipole-hadron 
Qh. (proton or nucleus) scattering amplitude at high rapidities [5, 6], different from the geometric scaling behaviour which 
is the hallmark of the "mean field" evolution equations (JIMWLK Q and BK [HQ equations), and which is named 
diffusive scaling by the authors of Ref. [l(| • The effect of fluctuations on the scattering amplitude [ll], [IH, [TH, [lj, [l5|, 
3 [H HHS, [2l|j2lj2|, the diffractive scattering processes [13, [H, SH and gluon production in hadronic 
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271 . 128j has been studied so far. In this work we focus on the consequences of fluctuations on 
the ratio R p a of the unintegrated gluon distribution of a nucleus over the unintegrated gluon distribution of a proton 
' scaled up by the atomic factor A 1 / 3 . 

The "Pomeron loop equations" describing the evolution of the dipole-hadron scattering amplitude with increasing 
rapidity Y = ln(l/x) are stochastic equations. In a frame where most of the rapidity is given to the hadron, the 
stochastic evolution of the hadron gives rise to different gluon distributions from one event to another 5j . The random 
variable in the evolution, the logarithm of the saturation momentum p s {A,Y) — hi(Q 2 (Y)/k 2 ), can therefore vary 
in an event by event basis. This variation is characterized by the dispersion a 2 = (p 2 ) — (p s ) 2 which rises linearly 
with rapidity, <J 2 (Y) = D& c a s Y, with D& c the dispersion coefficient and a s = a s N c /i:. For rapidities Y < Yds, the 
dispersion is small a 2 <C 1, meaning that the effect of fluctuations on the scaling form of the unintegrated gluon 
distribution can be neglected, while for Y > Yds, where a 2 ^> 1, fluctuations become important and do change the 
scaling form of the unintegrated gluon distribution. We will calculate the gluon distribution of the proton and the 
nucleus in the geometric scaling region at Y < Yds and in the diffusive scaling region at Y > Yds which are need in 
order to study R p a as a function of the atomic number A, the rapidity Y and the transverse gluon momentum k±. 
All calculations which are presented in this work are valid in the case of a fixed coupling a s . 
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In the overlapping geometric scaling regime of the proton and the nucleus RpA has been studied by using the 
JIMWLK or BK equations [H, M, M, HI \M, HI . It was found that R P A scales with A like A 1 / 3 ^ 0- ^ 

in the fixed-coupling case, rather than being equal to one (j c — 0.6275) as in the standard pQCD region at very 
large gluon momenta. This is partial gluon shadowing due to the anomalous behaviour of the unintegrated gluon 
distribution which stems from the BFKL evolution. Partial gluon shadowing may explain why particle production in 
heavy ion collisions scales, roughly like N part 37]. Furthermore R p a turns out basically k\ and Y independent in the 
common geometric scaling regime. The explanation is that the unintegrated gluon distributions of the nucleus and of 
the proton preserve their shapes with rising rapidity in the geometric scaling regime, yielding thus a constant value 
for their ratio, as shown for two different rapidities in Fig. HJa). We will derive the known results in the geometric 
scaling regime in this work in order to have a direct comparison with the results in the diffusive scaling regime. 

We have found that in the overlapping diffusive scaling regime of the proton and the nucleus R p a basically behaves, 
for a fixed coupling, as 
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where A.p s denotes the difference between the average saturation lines of the nucleus and the proton and (Q S (A, y)) 
is the average saturation momentum of the nucleus. This ratio shows two features which are different as compared 
to the ratio in the geometric scaling regime: (i) For k\ close to (Q S (A, T)) 2 , the gluon shadowing characterized by 

A 3 '^ -1 ' is dominated by fluctuations, through a 2 (Y), and depends also on the difference Ap s . Gluon shadowing 
increases as the rapidity increases because of a 2 = Dd c ct s Y. At asymptotic rapidity, where a 2 — > oo, one obtains 
total gluon shadowing, R p a = 1/A 1 ^ 3 , which means that the unintegrated gluon distribution of the nucleus and that 
of the proton become the same in the diffusive scaling regime. Total gluon shadowing is an effect of fluctuations since 
the fluctuations make the unintegrated gluon distributions of the nucleus and of the proton flatter and flatter [5[ and 
their ratio closer and closer to 1 (at fixed Ap s ) with rising rapidity as illustrated in FigfSJb). (Total gluon shadowing 
is not possible in the geometric scaling regime in the fixed-coupling case since the shapes of the gluon distributions 
of the nucleus and of the proton remain the same with increasing Y giving for their ratio a value unequal one, see 
FigEH a ))- (h) RpA shows an increase with rising k\ (always R p a < 1) within the diffusive scaling region. Since the 
exponent Ap s /a 2 decreases with rapidity, the slope of R p a as a function of k\ becomes smaller with increasing Y . 
The behaviour of R p a as a function of k± with increasing rapidity in the diffusive scaling regime is shown in Fig. [3] 

The JIMWLK [7] or BK [8, 9] equations, with the geometric scaling being their main consequence, appear to be 
appropriate for the understanding of the physics explored at RHIC and HERA experiments. However, it may be 
that this isn't the case anymore at LHC energies where fluctuations may start becoming nonnegligible. Therefore, 
in addition to the theoretically interesting results for R p a as a consequence of fluctuations, our result for R p a may 
also become relevant in the range of LHC energy. If this is the case, then the increase of the gluon shadowing and 
the decrease as a function of the gluon momentum of R p a with rising rapidity as given by Eq. ([T]) may be viewed as 
signatures for the onset of fluctuation effects in the LHC data. 

This work is organized as follows: In Sec.|H]we show two different definitions of the unintegrated gluon distributions 
used in the literature, define R p a and discuss the initial conditions for the proton and nucleus gluon distributions. In 
Sec. IIIII we show the results for the gluon distribution of the proton and the nucleus in the geometric and diffusive 
scaling regime. In Sec. IIVI we briefly review the known results for the ratio R p a in the geometric scaling regime, then 
calculate and discuss the new results for R p a in the diffusive scaling regime. 



II. GENERALITIES; DEFINITIONS, INITIAL CONDITIONS 



In this section we show two different definitions for the unintegrated gluon distribution used in the literature, define 
the ratio R p a and describe the initial conditions we use for the unintegrated gluon distribution of the proton and the 
nucleus. 



A. Definition of unintegrated gluon distributions and R p a 

Both definitions of the unintegrated gluon distributions of a hadron (proton or nucleus) known in the literature 
are expressed in terms of the forward scattering amplitude N (x±, bj_, Y) of a QCD dipole of transverse size x± with 
rapidity Y = In 1/x scattering off a hadron at impact parameter b± . Hereafter we consider the scattering process at 
a fixed impact parameter and omit therefore the fe^-dependence in all the following formulae. 
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The following definition has been proposed in Ref. [38| for the unintegrated gluon distribution: 

N f 

h A (k ± ,Y) - — I— / d 2 x ± e lk ^ Vl ± N(x ± ,Y), (2) 
(27r) a s J 

^fclvL / ^ e*-- TV (xx, y) • (3) 
7r a., J x 



(2 

This formula usually appears in cross sections for gluon production in proton- nucleus collisions [35j ] . The inversion of 
Eq. (0) being 

(27r) 3 a s /■ d 2 k ± , Ak,-x ± 



N ( X±7 Y) = ^ — - / h A (k ± ,Y) (2-e ik± - x± -e- ih ^) (4) 

V ' 2N C J (2nfkl 7 W 

does show more explicitly the relation between h A (k±,Y) and the dipole- nucleus scattering amplitude: The phase 
factors in the bracket are due to the four graphs describing the different ways of gluon exchange between the quark 
and antiquark of the dipole and the hadron. 

The other definition of the unintegrated gluon distribution reads 

( p A (k x ,Y) = -^l— [ ^Le^ k —N{x^Y) (5) 
(27t) a s J x ± 

and is derived from the non-Abelian Weizsacker- Williams field of a nucleus [H, H3, 5l| . This distribution counts the 
number of gluons in the wavefunction of the hadron. The two distributions are related to each other by 

h A (k ± ,Y) =k 2 ± V 2 kx <fi A (k ± ,Y) . (6) 

The main result of this work, the ratio R pA in the diffusive scaling regime, turns out to be basically the same for both 
distributions, as shown in Sec. lIVI Some more elaborate discussion on the two different definitions of the unintegrated 
gluon distribution can be found in Rcfs. j30l l33j| . 

The quantity which we study in this work is the ratio of the unintegrated gluon distribution of a nucleus over the 
unintegrated gluon distribution of a proton scaled up by A 1 / 3, , 



Rh h A (k ± ,Y) ^ <p A { k± ,Y) 

p Ai hp(k_L,Y) p Ai ip p (k±,Y) 



which is a measure of the ratio of the number of particles produces in a proton-nucleus collisions and the corresponding 
number in proton-proton collisions scaled up by the number of collisions [U HH, [H, US HH • 



B. Initial condition in the case of a nucleus 

The evolution of the unintegrated gluon distributions is given by the QCD evolution equations. However, one has 
to know the initial condition for the unintegrated gluon distribution to start with in the evolution equations. In the 
case of a nucleus, as an initial condition we use the McLerran-Venugopalan model [44| at fixed b± and Y = 0, 

N MV (x±, Y = 0) = 1 - exp L *± Q ^ A) ^j , (8) 

with the saturation momentum given by [44j 

Q 2 s {A) = 2 ^pT{b) X G{x,^j (9) 

where p is the nuclear density, T (b) = 2 yj R 2 A — b 2 is the profile function and xG (cc, is the gluon distribution 

in the nucleon which at Y — can be written as xG (x = 1,1/a;^) oc a s h\(\ / x^K^q D ) . Since we are interested 
throughout this work in the region where 1 / x± is not much larger than the saturation scale Q 2 (A) , we neglect the x± 
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by replacing it by 1/Q 2 S {A) in xG(x,l/x 2 L ). Noting also that pT(b) oc QqA 1 ^ 3 at fixed impact parameter, one gets 
for the a s and A dependence 



} 2 s (A)*Q 2 a 2 s Ai \n{a 2 A^) 



(10) 



which has also been used in Refs. [23,[3y|. We shall always assume that a 2 A^ ^> 1 in order to have non-trivial nuclear 
effects. 

Using the Mellin transform of the scattering amplitude, 



N M v (7) 



duu 7 1 [1 - exp (—it)] = — r (-7) with < 7 < 1 



where u = 



— — -, one can write Eq. ([8]) also in the form 

N M v (x±, Y = 0) = / p-N M v (7) exp 

J 27TI 



7 In 



QUA) 



(11) 



(12) 



which together with Eq. ([2]) and Eq. ([5]) allows us to get the unintegrated gluon distribution of a nucleus at Y = 0, 



h A (k±,Y = 0) 



•PA (k±,Y = 0) 



Nr. 



2n 2 a s \Q 2 (A) 



exp 



Q 2 S (A) 



Nr. 



8n 2 a 



■r 0, 



(13) 
(14) 



where F 



(°- A)) - ln (ofte) -7B 



is the incomplete Gamma function. 



Let us also construct the expressions for unintegrated gluon distributions which take into account BFKL evolution 
since they turn out to be useful in the next sections. This is easily done by starting with the BFKL evolved scattering 
amplitude in terms of the Mellin transform 



N M v (x±,Y) = 



c+zoo 



eh/ 



2ni 



-N M v (7) exp 



,Yx(l) +7 m 



AQUA) 



(15) 



where < c < 1, a s = as ^ a and x(l) — 2'0(1) — 4>(l) ~ ^(1 — 7) with ^(7) a polygamma function. With Eq.|T 
inserted into Eq.(f5]) and Eq.© and the relation 



d 2 *x Jk± . x± fxjQ 2 s (A) 



r(7) 



r(- 7 + i) \Q 2 M) 

one obtains the Mellin representation of the unintegrated gluon distributions which contain rapidity evolution: 



h A (k ± ,Y) - 
PA(k ± ,Y) = 



2ir 2 a 
N, 



c+ioo 



(hj 



r (7 + 1) exp 



s J c — ioc 

C+200 



8tt 2 a s 



s J c — zoo 



2tti 

dl r(7) 

2ni 7 



asYxh) - 7 In 



1(A) 



exp 



a s Yx(l) ~ 7 In 



QUA) 



(16) 

(17) 
(18) 



In the limit of Y = Eq. |(T5]). Eq. JT7I) and Eq. JTS]) reduce, of course, to Eq.((T2]l. Eq-JTSJ) and Eq.([T4]) respectively. 



C. Initial condition in the case of a proton 



We use two different ways to describe the initial gluon distribution of the proton: (i) In the main body of this 
work we view the proton as a single color dipole (see also Ref. (36|) and (ii) in Appendix [X] we view the proton as a 
composite object and use the McLerran-Venugopalan model to describe its initial conditions (see also Ref. [33]). We 
will show that both point of views exhibit equally well the effect of fluctuations on R p a which is the main focus of 
this work. 

The scattering amplitude for a dipole of size x±_ scattering off a dipole of size x' ± at relative rapidity Y reads [43j 



N(x' ± ,x x ,Y) 



2 12 
it a s Xj_ 



e+200 



1 



G?7 

7 2(1_ 7 )2 



exp 



a sX (A)F- 7 ln^ 



(19) 
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where < c < 1. When Y = the above expression reduces to the dipole-dipole scattering amplitude at the two 
gluon exchange level 

N(x' x ,x±,Y = 0) = 2-kc? s t\ fl + ln^ (20) 

with r < being the smaller of x± , x' ± and r> the larger of x± , Xj^ . 

Following the same steps as in the case of a nucleus, one obtains for the unintegrated gluon distributions of a proton 
which take into account BFKL evolution 

h P (k±,Y) = 7ra 2 x' 2 [ C+lC °±L ^_ Jf^ exp[s s yx(7)-7ln(fcixl 2 )] ; (21) 

2ir l a a / c _ joo 2tti (l-7) 2 r(l-7) 

= 8^7™'^ / c _ ioo 2™ f(l- 7 )»f(R exp[a s y X (7)-7ln(fc^J] • (22) 

The Eqs. JT7J), ([HI) and Eqs. (J2TJ , ([22j) have the form of the formulae one has started with in Ref. [H| and Ref. HQ 
to study the geometric scaling region and the diffusive scaling region. These equations tell us what the effect of the 
different inital conditions is according to Ref. [5, 6, 43]: Different initial conditions for the proton as compared with the 
ones for the nucleus do lead to different saturation momenta, Q s {p, Y) unequal Q S {A 1 Y). The saturation momentum 
does slightly depend also on the definition for the unintegrated gluon distribution, h Pj A (&_L , Y) and ip p ,A(k±, Y). In 
the next section we do show the exact expressions of the saturation momentum of the proton and the nucleus. 

III. UNINTEGRATED GLUON DISTRIBUTION OF A HADRON 

In this section we focus on the unintegrated gluon distribution of a highly evolved hadron (nucleus or proton) in 
the geometric and diffusive scaling region. To explain the relevant physics in these two regions let us look at the 
phase diagram of the hadron in the high-energy limit shown in the Y — p plane in Fig. [TJ Here Y = ln(l/a;) is 
the rapidity of the hadron and p = ln(fc^/feg) is the logarithm of the transverse momentum of gluons inside the 
hadron (k% is a fixed reference scale). The straight line denoted by (p s (A,y)) represents the average saturation line, 
(p s (A, Y)) = (ln(Q s (A, Y")/fc§)), with the saturation momentum Q 8 (A, y) depending on Y and the atomic number A. 
To the left of the saturation line, p < (p s (A,Y)), the gluon density is of order l/a s , or, the dipole-hadron scattering 
amplitude is at the unitarity limit (N) « 1. In this region, called "saturation region", the nonlinear QCD evolution 
becomes important. For p 3> (p s (A, Y)), the gluon density is low and the standard 'double-logarithmic approximation' 
is applicable. In the low density regime neither saturation nor fluctuation effects are important and the unintegrated 
gluon distribution behaves as l/fej. a ^ very large k± (the amplitude shows color transparency, (N(r±)) cx 7j_). The 
most interesting region for us in this work is the transition region between the saturation and the low density regime 
where p is not much larger than (p s (A, Y)) (see Fig. [TJ. There are two different re gim es within the transition region 
which are separated by the rapidity scale Yds, the geometric scaling regime pi |42|. [43[ and the diffusive scaling 
regime [||, in which the dynamics of the QCD evolution is different. 

The QCD evolution of the dipole-hadron scattering amplitude with increasing rapidity Y = ln(l/a;) is stochastic [H, 
6|] meaning that the scattering amplitude can fluctuate from event to event. Equivalently, in a frame where most of 
the rapidity is given to the hadron, the stochastic evolution of the hadron gives rise to different gluon distributions 
from one event to another [5[. Therefore the random variable, p s (A,Y) — \n(Ql(A, Y)/fcg), can vary from one event 
to another. Its variation is characterized by the dispersion a 2 = {p 2 s ) — (p s ) 2 which rises linearly with rapidity, 
<J 2 {Y) = Dd c a s Y , with D^c the dispersion coefficient. For rapidities Y < Yds, the dispersion is small a 2 <C 1, 
meaning that the effect of fluctuations on the scaling form of the unintegrated gluon distribution can be neglected, 
while for Y > Yds, where a 2 3> 1, fluctuations become important and do change the scaling form of the unintegrated 
gluon distribution as we show below (see also Ref. [l(| H3] for more discussions on the phase diagram) . 

A. Geometric scaling regime 

In the geometric scaling regime at Y <C Yp (Yp will be fixed below, see also Fig. [TJ the QCD evolution is described 
by the Kovchegov equation [9j , or, equivalently, by the BFKL equation in the presence of a saturation boundary 0, [43| . 
Starting with Eq. (fTTj) and following the calculations in Refs. 0, [43[ one finds for the unintegrated gluon distribution 
of a nucleus in the case of a fixed coupling 



In 



Ql{A,Y) 



1 

7c _ 



(23) 
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FIG. 1: The phase diagram of the unintegrated gluon distribution of a proton (or nucleus) in the kinematical plane Y— hn(kj_/k$). 
The high density region (saturation region) at small gluon momenta is on the left hand side of the saturation line (p s (A,Y)), 
the low density region at very high gluon momenta is one the right hand side of shadowed area. The intermediate region 
separated by Yds divides the geometric scaling regime from the diffusive scaling regime. 



which is valid in the geometric scaling regime 



Ql(A) <kl< Q 2 S (A) exp 



y/*x" (Tfc)a a Y 



(24) 



and scales as a function of the dimensionless variable k 2 /Q 2 (A, Y). The observation of a similar scaling behaviour at 
the HERA data for DIS at x < 0.01 and any momentum [46( supports the theoretical findings. In the above equation 
Q 2 S {A 1 Y) denotes the saturation momentum of the nucleus 

Ql (A, Y) = c h A Ql (A) ex V^ asY l (25) 
[2a s Y X "hc)}^ 

where ca is a constant and the value of the anomalous dimension, 7 C = 0.6275, which comes from the BFKL dynamics 
in the presence of saturation, is fixed by x'ilc) — x ^"^ ■ In Eq. (|23p ft, max = 0(l/a s ) is the maximum of tiA{k±, Y) 
which happens when k\ = Q 2 {A : Y). 

The unintegrated gluon distribution <p A (kj_, Y), defined in Eq. ([5]), has in the geometric scaling regime exactly the 
same form as JiA(k±, Y) given in Eq. ([23]) . The only difference, steaming from the different initial conditions at Y = 0, 
see eq. (1171) and Eq. (|18[) . is a slightly different constant cf A in the saturation momentum. 

To get the parametrical estimates for the rapidities Yp and Yds , it is convenient to express the above results in the 
geometric scaling regime in terms of the p — Ink^/k 2 variable, 

h A (k ± ,Y) = hT x lc e-^r-^ A > Y » (p-p s (A,Y) + ^J (26) 

with the region of validity 

0<p- P s(A,Y) < v /4 X "( 7c )a s T (27) 



and 



p s (A,Y)=)n^^ + ^a s Y- JL In [2a s Y X "( lc )] . (28) 
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Gluon number fluctuations (Pomeron loops) 0, Q do change the results (l26t and (|28|) which emerge from the 
Kovchegov equation as the rapidity increases. These results remain valid so long as the gluon occupancy in the 
nucleus, n(k±,Y) oc h A (k±, Y)/a s , is much larger than one in the geometric scaling regime. When n(k±,Y) « 0(1) 
(or the scattering amplitude N(r, Y) w a 2 s n(k ~ 1/r, F)) the discreteness of the gluons has to be taken into account, 
allowing the occupancy to go below one only by becoming zero [H, @] . The gluon occupancy becomes of order one at 
the largest k± within the geometric scaling regime, or h A (k±,Y) w a s , when p — p s m — \n(l/a 2 ). According to 

Y F oc — \n(l/a 2 s ) . (29) 



Eg. ([27)1 . this happens at the rapidity 



(p s (A, Y)) = In -J + ( — — 2/i ; ^ ) a s Y . (30) 



For Y F <Y fluctuations do change somewhat p s (A, Y) given in Eq. H25D to H, Ej 

1(A) (xhc) ^icx"{icY 

k'o V 7c 21n 2 (l/a2) 

Fluctuations, however, do not change the shape of the unintegrated gluon distribution in Eq. (|26[) in the region 
Yp < Y <C Yds because of the following reason: Starting with a unique initial gluon distribution at Y = 0, the 
stochastic evolution generates an ensemble of them at rapidity Y, where each of the individual gluon distributions has 
the shape given by ([26]) . but the individual gluon distributions may differ from each other by translation (p s (A,Y) is 
random). Based on the relation between QCD evolution and reaction-diffusion models in statistical physics [5(, the 
fluctuations are taken into account by averaging over all individual distributions, 



(h A (p - p s (A, Y))) = J d Ps h A (p - p s (A, Y)) P( Ps - ( Ps )) , (31) 
where the probability distribution of p s (A, Y) is argued to have a Gaussian form [45| . 

for p-p s (A,Y)«7 c V , (32) 



P (^) - fR — 9 exp 
V zna^ 



(Ps - (Ps)f 



2a 2 



= » ,ca v rcj (34) 



with the variance 

o 2 = (p s (A, Y) 2 ) - (p s (A, Y)) 2 = D dc a s Y (33) 
and the diffusion coefficient at asymptotic energies given by 0, Q 

*" 4 7cX"(7c) 
31og 3 (l/a 2 

Now, if the variance is small, a 2 <C 1, which is the case as long as the rapidity is much smaller than 

Yds * , (35) 

then the Gaussian distribution in Eq. (|32p is strongly peaked and therefore the averaged gluon distribution nearly 
preserves the form of an individual distribution, 

(MP ~ P>(A Y))) w h A (p - ( Ps (A, Y))) (36) 

and shows geometric scaling, at least in the window < p— p s (A, Y) < ln(l/a 2 )/7 c . So, at Y < Yp the unintegrated 
gluon distribution in the geometric scaling regime (see Fig. [I} is given by (|26p with the saturation momentum given by 
(|28p while at Yp < Y < Yds the unintegrated gluon distribution remains approximately the same while the saturation 
momentum changes to the one given in Eq. (|30|) . 

The only change in the case of a proton as compared to that of the nucleus in the geometric scaling regime is 
the saturation momentum, (p s (A,Y)) — ► (p s (p,Y)), which is done by replacing Q 2 (A) in the expressions for 
(p s (A,Y)) by Cp :V Q 2 (p), where c^ ,v is a constant and Q 2 (p)x' 2 — [a s -y/ln(l/ a s )] 2 ^ lc . The saturation momenta of 
the proton and the nucleus are different due to different initial conditions which we have explained in Sect. Ill Bl and 
Sect. Ill CI The difference between them is 

A P ^ {Ps (AY))-{p s ( P ,Y )} =1»||^ 



c c 



h ' v Ql x'l a 2 s A 1 / 3 ln(a 2 AV3 



In " A /* u "^ 7" ^ S ^ _ L (37) 
£ v [a. Vln(l/a,)] 2 ^ 
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with c a constant. A different Ap s is obtained when the proton is viewed as a composite object, instead of a single 
dipole, and is described by the McLerran-Vnugopalan model, as done in Appendix [A1 



B. Diffusive scaling regime 

The diffusive scaling regime sets in when the variance is large, a 2 3> 1, which happens at Y 3> Yds- I n the diffusive 
scaling regime the fluctuations are very important and the unintegrated gluon distribution changes therefore a lot 
from one event to another, leading to an average unintegrated gluon distribution (h A (p — p s (A, Y))) in Eq. (|31[) which 
is quite different as compared to the event-by-event gluon distribution — Ps{A, Y)) as we show below. 

We use for the event-by-event unintegrated gluon distributions 

h (n n (A Y] ,_jhT x exp[-(p s (A,Y)-p)} for p < p s (A,Y), , 

h A{ p - Ps (A, Y)j - < ^ ^ { _^ {p _ ps{A y))] for p > y) , (38) 



<p A (p-p a (A, Y)) 



. -max 



t^max gX p | 



[±(p s (A,Y)-p) + l] for p< Ps (A,Y), 
- 7c (p-p s (A,Y))} for p>p s (A,Y) 



(39) 



In the above equations we have used for p > p s (A, Y) the leading contributions valid in the geometric scaling regime 
(see Sect. IIII A[) while for p < pJA, Y) (saturation regime) the expressions which come from solving the Kovchegov 
equation (see Appendix [Bl and [9j) . There is a difference between the Ha and <pa distributions in the saturation 
region coming from their different definitions in {5} and ([5]). (We have used for simplicity the geometric scaling result 
also in the low density regime, see Fig. [TJ where both distributions behave the same way, like l/kj_, since the main 
result in the diffusive scaling regime is insensitive to the low density region.) 

It is now easy to show that the average unintegrated gluon distributions are roughly described by Gaussians, 



(h A ( P - p s (A,Y))) 
( V a(p - Ps(A,Y))) ~ 
and are valid in the diffusive scaling regime, 



i. max 

V2ncr 2 



l + 7c 



exp 



^max a 3 



2V2n[p-(Pe(A,Y))] 



■ exp 



(p-( Ps (AY))) 2 
2a 2 

\fi - (Ps{A,Y))] 2 



2a 2 



a < p- (p s {A,Y)) < 7c a 2 



(40) 
(41) 

(42) 



Note that in the diffusive scaling regime the unintegrated gluon distribution scales as a function of the dimcnsionlcss 
variable (p — (p s )) / a(Y) which is different as compared to the geometric scaling (p — (p s ))- The diffusive scaling 
extends up to very large values of gluon momenta k\ since the window p — (p s ) <C o 2 increases with rapidity, see 
Eg . (|33"]) . The same result as the one in Eq. ([40]) was obtained in Ref. (27[, although in a different way, and was shown 
to give the cross section for gluon production in the diffusive scaling regime. 

The geometric (diffusive) scaling window of a proton as compared to the geometric (diffusive) scaling window of 
a nucleus at high rapidity is shifted to lower momenta due to the smaller saturation momentum of the proton (see 
([37)) ). The proton and the nucleus have a common geometric scaling region for rapidities Y < Yds- The common 



diffusive scaling for the proton and the nucleus at Y > Y^§ sets in when the diffusive scaling window of the proton 



overlaps with the one of the nucleus, 



which determines 



Y pA or 



D dc a s Y > 



(Aps) 2 



r + A Ps 

for a > Ap s 
for a < Ap s 



(43) 



(44) 



1 The McLerran-Venugopalan model also gives the same expressions as a function of rapidity for p < p s (A, Y). 
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IV. THE RATIO R pA IN THE GEOMETRIC AND DIFFUSIVE SCALING REGIME 

In the geometric scaling regime the unintegrated gluon distributions Ha{p — (Ps)) and <Pa{p ~~ (Ps)) show the same 
behaviour (see Sect. IHI A| . From Eqs. (|26|) and (|37|) one finds for the ratio R v '£ defined in Eq. |7]) 



R h >X(k ± ,Y,A) 



pA 



{Qs(A,Y))< 



Ah [(Q S ( P ,Y)) 2 

1 [alHalA^)] 



In | 


^(Q S (A,Y))*) 








In 











Cj^ C X 



Ql 



where we have used h A lax 



AiO--yc) a 2 ln(l/a s ) 



//'""•' and ipj ax - 



„ h ,v> 



In 








In 


y(Q,(A,Y))*) 


+ Ap s + ± 





(45) 



ip' p """ and the average saturation momentum (Q s ) defined via 

(p s ) = \n(Q s ) 2 /kg. For gluon momenta not too far from (Q S (A, Y)) the ratio R^ scales with A like A 1 / 3 ^ lc ~ i \ 
rather than being equal to one. This is partial gluon shadowing due to the anomalous behaviour of the unintegrated 
gluon distribution given in Eq. (fT7|) which stems from the BFKL evolution. This result, already derived in Refs. [33TI36]]. 



may explain why particle production in heavy ion collisions scales, roughly, like N part [371 ]. Moreover R p ^ reduces to 
a k\ and Y independent expression for k\ much larger than (Q s (A,y)) 2 , or h\(k\/ (Q S (A : Y)) 2 ) 3> Ap s . This comes 
from the fact that the unintegrated gluon distribution of the nucleus and of the proton preserves the shape with rising 
rapidity, yielding therefore a constant value for their ratio. This behaviour is shown in Fig. [2ja). 

In the diffusive scaling regime the ratio R p a for the unintegrated gluon distribution in Eq. (|40[) equals 



RpA(k±-, 



Y,A) 



1 


UQ S {AX)?~ 


2a- 2 


\ ^ 1 


II 


{{Qs(p,Y)) 2 _ 




[(Q S (A 7 Y)} 2 \ 



1 


' a 2 ln^A 1 / 3 ) " 


2cr' 2 




2a 2 


\ ^ 1 


Az (L ^ T> 


[(a 2 \n(l/a s )y/^_ 




L c p 




[(Q S (A,Y)) 2 \ 



(46) 



while for the " Weizsacker- Williams" unintegrated gluon distribution in Eq. (|4"Tj) it is 

2 



n pA 



1 


\{Q S (A,Y)) 2 ~ 


2a-' 2 


\ ^ 1 


H 


[(Q s (p,Y)) 2 _ 




[(Q S (A,Y)) 2 \ 



hi 



((Q 3 (A,Y)y) +A P* 



111 



\(Q S (A,Y)) 2 ) 



1 


' a 2 s ln(«2 AV3) " 


Ap s 
2a' 2 


C A c x ± Qo 


2 a- 2 


\ 1 


A :i[ ^ T > 


i(a2ln(l/ as ))i/7c_ 








[(Q S (A,Y)) 2 \ 



111 



(Qs(A,Y)y 



Ap s 



In 



(47) 



In the diffusive scaling regime R p A becomes independent of the definition of the gluon distribution and shows the 

universal behaviour v4 1/3( ^ _1) [k 2 _/(Q s (A, Y)) 2 ] when \n(k\/ {Q S {A, Y)) 2 ) > Ap s . The above ratios in the 
diffusive scaling regime show two different features as compared to the ratio in the geometric scaling regime: (i) For 

k\ close to (Q S (A,Y)) 2 , the gluon shadowing characterized by A 1 ^ 3 ^ ^ is dominated by fluctuations, through 
<7 2 (Y), and depends also on the difference Ap s . The gluon shadowing increases as the rapidity increases because of 
a 2 = Dd c a s Y. At asymptotic rapidity, where a 2 — > oo, one obtains total gluon shadowing, R'a — 1/A 1 / 3 , which 
means that the unintegrated gluon distribution of the nucleus and that of the proton become the same in the diffusive 
scaling regime. The total gluon shadowing is an effect of fluctuations at fixed couping since the fluctuations make 
the unintegrated gluon distributions of the nucleus and of the proton flatter and flatter and their ratio closer and 
closer to 1 (at fixed Ap s ) with rising rapidity, as shown in Fig. [2Kb) . (ii) R p A shows an increase with rising k\ within 

the diffusive scaling region. Since the exponent Ap s /a 2 decreases with rapidity, the slope of R p A as a function of 
k\ becomes smaller with increasing Y. Both features of R p a in the diffusive scaling regime, the increasing gluon 
shadowing and the descreasing k± dependence with rising rapidity, are shown in Fig. [3] Note also that in the common 
diffusive scaling regime, see (|42l) and (|4"3")) . the ratio R p a is always smaller than one. (This can be seen especially from 
Eq. (|A5[) where Ap s = In A 1 / 3, was used.) 

Total gluon shadowing is not possible in the geometric scaling regime in the fixed-coupling case since the shapes of 
the gluon distributions of the nucleus and of the proton remain the same with increasing Y giving for their ratio a 
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value unequal one (see Fig[2ja)). In the case of a running coupling, the gluon shadowing increases with rising rapidity 
in the geometric scaling regime [33[ , as oppossed to the (roughly) fixed value (partial shadowing) in the fixed-coupling 
case, and would lead to total gluon shadowing at very high rapidities if fluctuations were absent. The combination of 
the running of the coupling plus fluctuations would give total gluon shadowing at lower rapidities as compared to the 
fixed coupling case presented here. An extension of this work by the running coupling remains work for the future. 





FIG. 2: The qualitative behaviour of the unintegrated gluon distribution of a nucleus (A) and a proton (p) at two different 
rapidities in the geometric scaling regime (a) and diffusive scaling regime (b). 




FIG. 3: The qualitative behaviour of the ratio R p a as a function of k± at four different rapidities, Y\ < Y2 < Y3 < Y4, in the 
diffusive scaling regime. R p a is always smaller than one for values of k± in the diffusive scaling regime. 
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APPENDIX A: MCLERRAN-VENUGOPALAN INITIAL CONDITION FOR THE PROTON 

Also the proton can be described within the McLerran-Venugopalan model, as done for the nucleus in Sect. Ill Bl 
The only change as compared to the description of the nucleus is a different saturation momentum for the proton due 
to the different initial conditions. It is reasonable to assume the following relation between their saturation momenta 



{Q s (AX)f = A 1 '* (Q s (p,Y)) 2 



(Al) 
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or in terms of (p s (A,Y)) = ln((Q s (A, Y)) 2 /k 2 ), 

( Ps {A,Y)) = {p s {p,Y))+\nA 1 '* 



(A2) 



This relation holds formally if one extrapolates the McLerran-Venugopalan model, which was constructed for a large 
nucleus, down to A = 1. So, to get the expressions for the proton from the ones given for the nucleus in the previous 
sections, one has to replace (Q S (A, Y)) 2 by (Q s (p,Y)) 2 , or, the difference Ap s in Eq. ([37]) by 



Ap™ v = ( Ps (A,Y))-(p s (p,Y))=lnAi . 
With this difference, the ratio in the geometric scaling region, as compared to (jl5|) . now reads 



R h p 'X(k ± ,Y,A) 



1 



In 



Ai 



*(i-r„) ]u ( fcj \ 

111 {(Q„{A,Y))* ) 



+ In A s 



(A3) 



(A4) 



and shows the same features (partial gluon shadowing and k\ and Y independence for In k 2 _/(Q s (A, Y)) 3> ln^ 1 / 3 ) 
as the result in Eq. (|45|) . In the diffusive scaling regime, instead of Eq. (|46|) . now one obtains 



rth 



A T 



k\ 



{Qs{A,Y)f 



(A5) 



and instead of Eq. ([17)) 



H pA 



A" 



{Qs(AY)Y 



In A 1 / 3 

„2 



In ( 



(Qs(A,Y)V 



In An 



In 



( *3 ^ 



(A6) 



which show the same universal features for the ratio R'^ , total shadowing at 1" — > oo and a decreasing k\ dependence 
with rising Y, as discussed in the previous section. 



APPENDIX B: GLUON DISTRIBUTION IN THE SATURATION REGION 



We calculate the gluon distribution in the saturation region, where k\ < Q 2 S , by solving the Kovchegov equation 9] 
at high rapidities. The Kovchegov equation in momentum space, at a fixed impact parameter, reads [9( 



d<f>(p,Y) 
a x dY 



X {~d p )4>{p,Y)-cp 2 (p,Y) 



(Bl) 



where p = In ^ and (f>(p,Y) = ( 27 a » ^ A (k±,Y) the Weizsacker- Williams gluon distribution defined in Eq. §5$) 
(apart from the prefactor) . It is convenient to look for a solution to the Kovchegov equation in the saturation regime 
in the form 



4>(P,Y) 



c+i °° d-y 

—<j> (7) exp [-7 (p - va s Y)} 



c— ZOO 
c+ioo 



-^-A> (7) exp 



-7 In 



(B2) 
(B3) 



with Q 2 (Y) = k 2 exp [a s »y]. In the saturation regime, In ( c^fy) ) ^ ^' ^ ne con t° ur integral has to be closed to the 

left and form a counterclockwise path in order to enclose all the poles of <f> (7). The first and most important pole of 

00 

4> (7) is at 7 = 0. Using the Laurent expansion of <j) (7) at 7 = 0, which is <f> (7) = a nJ n , in Eq. ljBlj) . one can 

n— — 00 

immediately spot that there must be a cutoff of n at n — —2. Thus, <f> (7) can be written as 



0(7) = ^ + - + «3 + /(7) 

7 7 



(B4) 
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with /(O) = 0. The function / (7) would have some poles at 7 = — 7^ where 7$ > 0, which however do generate small 
contributions like (J^)" to <P(p,Y) in the saturation region. Therefore the most important contributions are 
from O (7) = <t>{l) - /(7)< 



It is now straightforward to see that 

^o(p,r) = 01 ln^|P 



a 2 , 



X (-dp) 0o Go, F) = - a x [ In -33 



Ql (Y) 



a 2 In 



(Y) 



a 3 



Eft-- 



n=l 



«3 



Qs (Y) 



(B5) 
(B6) 
(B7) 



where the last term in (|B7I) comes from the residues of x (7) at 7 = —n with n — 1,2,. 
correction of order „f jy. , Eq. ljBip is solved by O (p, Y) = a\ In ^ 



Qj(y) 



a 2 with ai = 5 and \v 



Up to some small 
= 03 — a\. Moreover 



one can set a 2 = and get 03 = ^v. Finally, one obtains the approximate solution to the Kovchegov equation, 
(p, Y) = I In + O ( q^(V) ) ' anc ^ ^he g mon distribution as defined in Eq. ([5|) takes the form in the saturation 



regime as cited in Eq.(|39 



<PA (k±,Y) = 



In- 



(9 



2(27r)^a s L *i W 
The gluon distribution as defined in Eq.@ reads in the saturation regime 



h A [k^,Y) oc 



2(2nfa s 0*00' 

and is obtained from Eq. (|B8|) , the higher order correction term, by using the relation in Eq. 



(B8) 



(B9) 
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